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Abstract 

To increase the performance of Stirling engines and analyze their operations, a second-order Stirling model, which includes thermal 
losses, has been developed and used to optimize the performance and design parameters of the engine. This model has been tested using 
the experimental data obtained from the General Motor GPU-3 Stirling engine prototype. The model has also been used to investigate 
the effect of the geometrical and physical parameters on Stirling engine performance and to determine the optimal parameters for 
acceptable operational gas pressure. When the optimal design parameters are introduced in the model, the engine efficiency increases 
from 39% to 51%; the engine power is enhanced by approximately 20%, whereas the engine average pressure increases slightly. 

© 2008 Elsevier Ltd. All rights reserved. 
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1. Introduction 

An elementary Stirling engine is composed of an engine 
piston, an exchanger piston, and three heat exchangers: a 
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Nomenclature 

CO 

operating frequency (rd/s) 

A 

area (m 2 ) 

Subscripts 

C p 

c 

^pr 

specific heat at constant pressure (J/(kgK)) 
heat capacity of each cell matrix (W/K) 

c 

compression space 

£ 

regenerator efficiency 

d 

expansion space 

h 

convection heat transfer coefficient (W/(m 2 K)) 

diss 

dissipation 

M 

mass of working gas in the engine (kg) 

E 

entered 

m 

m 

mass flow rate (kg/s) 

mass of gas in different components (kg) 

f 

cooler 

P 

pressure (Pa) 

h 

heater 

Q 

heat (J) 

irr 

irreversible 

Q 

power (W) 

P 

loss 

R 

gas constant (J/(kgK)) 

Pa 

wall 

T 

temperature (K) 

r 

regenerator 

k 

wall conductivity (W/(mK)) 

S 

outlet 

V 

volume (m 3 ) 

shtl 

shuttle 

w 

work (J) 

T 

total 


as shown in Fig. 1. The engine uses external combustion, 
hence it can be powered by any source of energy 
(combustion energy, solar energy, etc.) and causes less 
pollution than the traditional engines [1-5]. The working 
piston converts gas pressure into mechanical power, 
whereas the exchanger piston is used to move gas between 
hot and cold working spaces. The engine presents an 
excellent theoretical efficiency of the same order as the 
Carnot efficiency. 

Several prototypes have already been produced [2-5], 
but their actual efficiency remains very low compared to 



the high theoretical efficiency. In fact, Stirling engines 
involve extremely complex phenomena related to compres¬ 
sible fluid mechanics, thermodynamics, and heat transfer. 
An accurate description and understanding of these highly 
non-stationary phenomena is necessary so that different 
engine losses and optimal design parameters may be 
determined. 

Several authors have studied the finite-time thermody¬ 
namic performance of Stirling engines and the effect of 
heat losses and irreversibilities on engine performance 
[6-11]. However, they have not calculated the optimal 
design parameters for maximum power and efficiency. 
Popescu et al. [6] have shown that the low performance is 
mainly due to the non-adiabatic regenerator. Kaushic, Wu, 
and co-workers [7,13,15] have proved that the most 
important factors affecting the performance of Stirling 
engines are heat conductance between the engine and 
reservoirs, the imperfect regenerator coefficient and the 
rates of the regenerating process. Kongtragool and 
Wongwises [9] have investigated the effect of regenerator 
efficiency and dead volume on the engine network, heat 
input, and engine efficiency, using a theoretical approach 
to the thermodynamic analysis of Stirling engines. Costea 
et al. [10] have studied the effect of irreversibility on solar 
Stirling engine cycle performance. They have included the 
effect of incomplete heat regeneration, and internal and 
external irreversibility of the cycle. Cinar et al. [12] 
manufactured a beta-type Stirling engine operating at 
atmospheric pressure. The tests carried out on this engine 
have shown that the engine speed, engine torque, and 
power output increase with the hot source temperature. 
Walker [2] has identified several other losses, such as 
conduction losses in the exchangers, dissipation by pressure 
drop, shuttle and gas spring hysteresis losses. These losses, 
however, are not usually accounted for in the published 
work due to their complexity. Urieli and Berchowitz [14] 
have developed a quasi-steady flow model that includes 
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only the pressure drop in the exchangers. The results 
obtained by the author using this model are more accurate 
than those obtained by other models, though they are 
different from corresponding experimental data. Abdullah 
et al. [16] have given the design considerations to be taken 
into account when designing a low-temperature differential 
double-acting Stirling engine for solar applications. They 
have determined the optimal design configuration, such as 
the engine speed, regenerator porosity, and heat exchanger 
volumes, including the pressure drop, heat exchanger 
effectiveness, and the swept volume. Therefore, the Stirling 
engine performance depends on geometrical and physical 
parameters of the engine and the working fluid, such 
as regenerator efficiency, porosity, dead volume, swept 
volume, temperature of sources, energy and shuttle 
losses, etc. 

A numerical simulation that accounts for losses has been 
developed by the authors and tested using the General 
Motor GPU-3 Stirling engine data [17]. The results 
obtained proved better than those obtained by other 
models and correlates more closely with the corresponding 
experimental data. The model is used to determine losses in 
different engine compartments and to calculate the 
geometrical and physical parameters corresponding to 
minimal losses [18,19]. An optimization based on this 
model is presented in this article. This study helps to 
determine the influence of geometrical and physical 
parameters on prototype performance and therefore to 
identify the optimal design parameters. 

2. Dynamic model with losses 

A second-order adiabatic model has been initially devel¬ 
oped; the estimated engine parameters are computed for use 
in the validation. The results are compared with those 
obtained by Urieli and Berchowitz [14] at the same conditions 
(adiabatic model). Afterward, a dynamic model including 
losses in different engine elements has been developed. 

The losses considered in this model include energy 
dissipation by pressure drops and internal conduction 
through exchangers. Furthermore, the following energy 
losses have also been included: the external conduction in 
the regenerator, the shuttle effect in the displacer, and the 
gas spring hysteresis in the compression and expansion 
spaces. However, mechanical friction between moving 
parts has not been considered. 

2.1. Losses included in the model 


The losses considered in the model and evaluated in 
articles [17,18], are: 


1. energy dissipation by pressure drops in heat exchangers 
(SSdissX given by [14] 



(1) 


where Ap, the frictional drag force, is given by 


Ap = 


2 fdiGV 
Ad 2 p 


( 2 ) 


where G the working gas mass flow (kgm~ s~ ), d the 
hydraulic diameter, p the gas density (kgm - ), V the 
volume (m ), and / r the Reynolds friction factor; 

2. energy loss due to internal conduction (8(2 pc d) between 
the hot parts and the cold ones of the engine through 
the different exchangers [14]: 


§&> 


c d rl 


G d rl (2"r—r “ G-r) 


(3) 


8(2pcdr2 — Gdr2 j (G-h G-r) 


(4) 


§(2pcdf = Gdf y~(G-r - G-f) 
G 


(5) 


^0Pcdh — Gdh j (G-d G-h) 
G 


( 6 ) 


where k cd (W m -1 K -1 ) is the thermal conductance of 
the material and A the effective area for conduction; 

3. energy loss due to external conduction (8£> pe xt) in the 
regenerator which is not adiabatic. These losses are 
specified by the regenerator adiabatic coefficient s^l, 
defined as the ratio between the heat given up in the 
regenerator by the working gas during its transition 
toward the compression space and the heat received in the 
regenerator by the working gas during its transition 
toward the expansion space [3]. Hence the energy stored 
by the regenerator at the time of the transition of the gas 
from the expansion space to the compression space is not 
completely restored with this gas at the time of its return. 
For the ideal case of the regenerator with perfect in¬ 
sulation, s = 1. The energy lost by external conduction is 


Sgpext = (1 - eXS&i + § Q i2 ) 


(7) 


The effectiveness of the regenerator ds given starting 
from the equation below [8,14] : 


NTU 
1 +NTU 


( 8 ) 


NTU is the number of transfer units: 
NTU — g 


C p m 


(9) 


where h is the overall heat transfer coefficient (hot stream/ 
matrix/cold stream), A wg refers to the wall/gas, or 
“wetted” area of the heat exchanger surface, C p the 
specific heat capacity at constant pressure, and rh the mass 
flow rate through the regenerator. 

4. energy loss due to shuttle effect (8£) psh ti): the displacer 
absorbs a quantity of heat from the hot source and 
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restores it to the cold one. This loss of energy is given 
by [16] 

?ih _ 0.4Z 2 /q,i s T> d T . /i r\\ 

8£3p s htl — Od-lc) (lb) 

where J is the annular gap between the displacer and the 
cylinder (m), k pis the piston thermal conductivity 
(Wm _1 K _1 ), D d the displacer diameter (m), T d the 
displacer length (m), Z the displacer stroke (m), and T d 
and T c are, respectively, the temperature in the 
expansion space and in the compression space (K); 

5. energy loss due to gas spring hysteresis (5(7 P i rr ): For an 
ideal gas, the pressure/volume relationship is either 
isothermal or adiabatic. In a real gas, there is a certain 
amount of work that is dissipated (S£)pin- = —5fF irr ). 
Urieli and Berchowitz [14] gave the expression of this 
loss. For the compression space we have 


iVFh 

Rmu 


(15) 


PdVd 

Rm d 


(16) 


(4) the regenerator is divided into two cells rl and r2; each 
cell has been associated with its respective mixed mean 
gas temperature T rl and T t2 expressed as follows: 


Pn Frl 

Rm^\ 


(17) 


Pr2 F r2 
Rm V 2 


(18) 


1 


6 IFp j rr c — ^ COy^iy l)Tp ac T cm0 y/t Pac | 


jk\ 


A k r 


V, 


A p. 


c moy y 


Pa c 


(ii) 


For the expansion space we have 


5 W 


P irr d 



l)Tp a d P d moykpa 


’ A Vj V 

Vi moy / 


^4 Pa d 

( 12 ) 


where y — C p /C v , T Pac and T Pad are the wall tempera¬ 
tures in the compression space and the expansion space, 
respectively, which are equal to the average temperature 
in these spaces; AV C and AF d are, respectively, the 
volume amplitude in the compression space and the 
expansion space, ^4 Pac and ^4p ad are the wetted area in 
the compression space and the expansion space, 
respectively, and co is the operating frequency. 


2.2. Model development 

The schematic model of the engine and various 
temperature distributions in the engine components are 
shown in Fig. 2. The dynamic model of the developed 
Stirling engine is based on the following assumptions: 

(1) the gas temperature in the different engine elements 
varies linearly; 

(2) the cooler and the heater walls are maintained 
isothermally at temperatures T Paf and T Pah ; 

(3) the gas temperature in the different components is 
calculated using the perfect gas law: 


An extrapolated linear curve is drawn through tempera¬ 
ture values T r \ and T r2 , defining the regenerator interface 
temperature T f _ r , T r _ r , and T r _ h , as follows [18]: 


T f- r 


3T r i — T r 2 
2 


(19) 


T 

J. r—r 


T rl + Tv2 
2 


( 20 ) 


Tr -h = 


3T r 2 — T r i 
2 


( 21 ) 


According to the flow direction of the fluid, the 
interface’s temperatures shown in Fig. 2 are defined as 
follows [17]: 

T c _ f , the temperature of the interface between the 
compression space and the cooler, is T c _ f = T c if m c _ f >0, 
otherwise T c _ f = T f . 

7f_ r , the temperature of the interface between the cooler 
and the regenerator, is 7f_ r = T { if m f _ r >0, otherwise 

Tf-r = Ty- f. 

T r _ h , the temperature of the interface between the 
regenerator and the heater, is T r _ h = 7i-_ r if m r _ h > 0, 
otherwise T r _ h = T h . 

7h_ d , the temperature of the interface between the heater 
and the expansion space, is 7h_ d = Th if m h _ d > 0, 
otherwise Th_ d = T d . 

Heat is transferred to the working gas by means of 
forced convection given by 


§2 = /^Pa(T Pa - Tf) 


( 22 ) 


P C V C 

Rm c 


Rmf 


where h is the heat transfer coefficient. Taking into account 
the losses by internal conduction in the exchangers and 
external conduction in the regenerator, the power ex¬ 
changed in the different heat exchangers are given by 

(14) 8 Q t = MpafUpaf - Tf) - 5e Pcdf 


(23) 
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Compression 

space 



m c-f 

ivT 


cooler 

Vf 
» m f 

Tf 


Tpaf 


Of 


Regenerator 


Expansion 

space 


"rl ,m,i V f 2 ,m r 2 
Ttl.Prl Ti2. p f2 


“f.r %-r m r-h 

Tjp. f " Tf-r T f -h 


Heater 


In h 

Th 


- / 

VQrl 

/ 

rO 

to 

Hi 





T Pail Tp^ 


A 


Qh 


f'h 

T Pah 


W d 


% 

m h-d m c j 

T h-d T d 

— Pd 



S&, = e/ lrl ^ Pad (r Parl - J rl ) - Sgp 2 cdrl (24) 

se r2 = sh r 2 Ap a r2 (r Pa r2 - r r2 ) - sgp 2 cdr2 (25) 

§2h = ^4pah(7pah — E h ) “ §2p c dh (26) 

where c is the regenerator efficiency. 

The heat transfer coefficient of exchanges, hf, h ri , h r2 , 
and Ah, are available only empirically, being complicated 
functions of the fluid transport properties, the flow regime, 
and the heat exchanger geometry [14]. 

The regenerator matrix temperatures are, therefore, 


given by 



dTp ar i 

52rl 

(27) 

dt 

Cp r dt 

d T Pa r2 

52 r2 

(28) 

dt 

Cp r dt 


There is no leakage, the total mass of gas in the system 
( M) being constant. Thus 


M — m c + mf + m r \ + m r 2 + (29) 

The energy equation applied to a generalized cell is 
reproduced as follows: 

§e + C„ T E m E - Cp r s »is = 4 W + C v (30) 


Applying the energy equation to the compression space, 
we obtain 


50c - C p r c _ f w c _f = ^ W c + (31) 

Since the compression space is adiabatic, 52c = 0 and the 
work done is dfE c /dt = iME c /dt. From continuity con¬ 
siderations, the rate of accumulation of gas (w c ) is equal 
to the mass inflow of gas, given by —m c _f. Thus, Eq. (31) 
reduces to 


C p T c _f m c = P ^ E c + C v ^ (m c T c ) (32) 

Substituting the equation of state and the associated ideal 
gas relations C p —Cy = R, C p = Ry/y— 1, and Cy = 1 

into Eq. (32) and simplifying gives 


th c 


1 (f dVc | AdP\ 

«r c _ f V d/ j di) 


(33) 


Taking into account the loss of gas spring hysteresis 
in the compression and expansion space, 5JEi rrc /dt and 
8fEj r rd/cU, the work generated by the cycle can be expre¬ 
ssed as 


8fV_ dV c dV d 5 W irrc 5 FE irr d 

dt c dt d dt d t d t 

The total engine volume is 

E t = E c + E f + V r \ + E r2 + E h + E d 


(35) 
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Since there is a variable pressure distribution throughout 
the engine, we have arbitrarily chosen the compression 
space pressure P c as the baseline pressure. For each 
increment of the solution, P c is evaluated from the relevant 
differential equation and the pressure distribution is 
determined with respect to P c . Thus, it can be obtained 


from the following expression: 


AE f 

P[ = p c+ 2 f 

(36) 

(APf + Aiy) 

P, i -P f + 2 

(37) 

D D , (AE r i E AE r2 ) 

G 2 = Gi E ^ 

(38) 

P^Prl + ^’Q^ 

(39) 

P.-P. + ^ 

(40) 


The other variables of the dynamic model with losses are 
given by the energy and mass conservation equation 
applied to a generalized cell. Taking into account energy 
dissipation caused by pressure drop in the exchangers 
(§£?diss) and the other losses yields 

dQ - 5<2 diss - b^Pshti + C p T E m E - CpTsths 

dW d W irr , ^ d (mT) (41) 

d t dr v dr 


Substituting the ideal gas relations into Eq. (41) and 
simplifying gives 


bQ - bQ dlss - bQ PsM + C p T E m E - C p T s ms 


1 

R 



E CyE 



-bW 


irr 


(42) 


Applying expanded energy conservation Eq. (42) to the 
different engine cells in Fig. 2 gives 


CpT c —{fh c —{ — 

K 


CvPcG+CvKcG) -8 W 


dt 


dt 


irr c 


(43) 


§(2f — §2diss f + CpTVfWc-f — CpTf_ r mf_ r = 


§(2rl ^2dissrl "E CpT f_ r /hf_ r Cp T r _ r /U r _ r — 


dQr2 ^(2dissr2 "E CyT Cp T r — E m r — E — 


Cy Vfdp c 

R 

dt 

(44) 

Cy Erl 

dP c 

R 

dr 

(45) 

Cy V r 2 dP c 

R 

dr 


( 46 ) 


dQ E ^2dissh E CpTf— E TVl r —h Cp Th—dfilh—d 


CyV h dP c 
R dt 
(47) 


1 


dV d 


CpTh—d^h—d bQ P sht i — I (TpP d ECv Ed 


R 


dr 


dA 

dr 


- bW 


irr d 


(48) 

Summing Eqs. (43)-(48), we obtain the pressure variation: 


d P c 


1 


dr CyV E 


R(bQ - bQ dissT ) - C l 


bW' 

dr 


(49) 


where bQ = 5(9f+ bQ rl + bQ r2 + bQ h -bQp shtl is the total 
heat exchanged. 6£7dissT 6£7dissfE §2dissrl E §2dissr2E 
§£?dissh is the total energy dissipation generated by pressure 
drop. 

There is no flow dissipation in the compression space; 
the mass flow of gas, Eq. (33), remains unchanged: 


m c = 


_J_/p 

RT c -f V c dr c 7 dr 


(50) 


The mass flow in the different engine components is 
given by the expanded energy conservation equations 
(43)-(48): 


1 


m c - f = - 


m-r = 


RT c - { 
1 


dEc 

dr 


+ V c 


dPQ 
ydt , 


+ 


6 IV 


irr 


C n G- 


(51) 


p 1 c-f 


m r— r = 


m T - h = 


m h - d 


CpT{_ r ' 

-J-( 

Cp T r _y \ 

Cp T r _h 

1 

C p Th_d 


bQf 5(2diss f + Cp T c _fm c _f 


bQr 1 ^Qdissrl "E Cp Tf_ r /hf_ r 


bQ r2 5Qdissr2 E Cp T r _ r fh r _ r 


bQy ^Qdiss h E C p T r _ h m r _ 


Cy Vf dP c \ 

R dt J 

(52) 

Cy V r \ d P c \ 
R dt J 

(53) 

Cy V r 2 dP c \ 
R dt J 

(54) 

Cv F h dP c \ 
R dt J 


(55) 

The equation of continuity is recalled as follows: 
m = m E — riis (56) 

Successively applying Eq. (56) to the four heat exchanger 
cells in Fig. 2, we obtain 

thf — rh c -f — thf-r 


m r 1 = Wf_ r — m r _ r 
m r2 = m r -r - w r _h 
m h = m r _ h - m h _ d 


(57) 

(58) 

(59) 

(60) 


3. Method of solution 

The model developed has been tested using data from the 
Stirling engine GPU-3 manufactured by General Motor. 
This engine has a rhombic motion transmission system, as 
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shown in Fig. 1. The geometrical parameters of this engine 
are given in Table 1. The operating conditions are 
as follows: working gas helium at a mean pressure of 
4.13 MPa; frequency 41.72 Hz; hot space temperature 
T Pah = 977K; cold space temperature T Paf = 288K. The 
measured power output was 3958 W, at a thermal efficiency 
of 35%. 

The independent differential equations obtained in 
paragraph 2, are solved simultaneously for the variables 
P c , m c , T rU W , etc. The vector Y denotes the unknown 
functions. For example, Y Pc is the system gas pressure in 
the compression space. The initial conditions to be satisfied 
are noted: Y(t 0 ) = F 0 . 

The corresponding set of differential equations is 
expressed as dY/dt = F(t, Y). The objective is to find the 
unknown function Y{t) which satisfies both the differential 
equations and the initial conditions. The numerical 
solution is composed of a series of short straight-line 
segments that approximate the true curve Y{t). It starts 
from the stationary state, with T c and T d at F Paf and F Pah 
or any arbitrary initial temperature values, and goes 
through successive transient cycles until the values of all 
the state variables at the end of each cycle are equal to their 
values at the beginning of that cycle. The system of 
equations is solved numerically using the classical fourth- 
order Runge-Kutta method, cycle after cycle until periodic 
conditions are reached. 

To validate the numerical method used in the computa¬ 
tion, the results are compared with those obtained by Urieli 
and Berchowitz [14] for the same conditions (adiabatic 
models) of the GPU-3 engine data. The comparison shows 
a good agreement, as shown in Fig. 3. 

4. Dynamic model results 

The model has been developed gradually, initially by 
neglecting the losses, then by introducing them progres¬ 


sively. The comparison of the various models results with 
those obtained by Urieli and Berchowitz [14] for the same 
conditions of the GPU-3 engine data is shown in Table 2. 
The comparison shows a good agreement. 

To show the effect of each loss on the engine’s 
performances, we represented the results of the model by 
gradually integrating the various losses. 

When all losses (50 Pcdf , §2dissh, §<2p sh ti, etc.) are 
included in the model, the heat flow rate for each compo¬ 
nent versus crank angle is illustrated in Fig. 4. The corres¬ 
ponding average power of the engine is equal to 4.27 kW. 
The average heat flow generated by the heater is equal 
to 10.8kW; it leads to an engine efficiency of 39.5%. 



- B erchowitz re suits (Adiabatic model) 

X X X Our results (Adiabatic model) 

Fig. 3. Validation of the computational method. 


Table 1 

Geometric parameter values of the GPU-3 Stirling engine 


Parameters 

Values 

Parameters 

Values 

Clearance volumes 


Cooler 


Compression space 

28.68 cm 3 

Tube number/cylinder 

312 

Expansion space 

30.52 cm 3 

Interns tube 

1.08nnn 

Swept volumes 


Diameter 

46.1 mm 

Compression space 

113.14 cm 3 

Length of the tube 

13.8 cm 3 

Expansion space 

120.82 cm 3 

Void volume 


Exchanger piston conductivity 

15 W/mK 

Regenerator 


Exchanger piston stroke 

46 mm 

Diameter 

22.6nnn 



Length 

22.6nnn 

Heater 


Wire diameter 

40 pm 

Tube number 

40 

Porosity 

0.697 

Tube inside diameter 

3.02nnn 

Unit number/cylinder 

8 

Tube length 

245.3 mm 

Thermal conductivity 

15 W/mK 

Void volume 

70.88 cm 3 

Void volume 

50.55 cm 3 
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The power and the efficiency calculated by the model are 
very close to the power and the actual efficiency of the 
prototype given in the abstract, but compared to those of 
Martini [20], we note that they are a little different from 
those of Martini tests 5 and 6. This is probably due to 
different operating conditions and setting of equations. 

The heat flow lost by internal conduction, the energy 
dissipation by pressure drop through the heat exchangers 
and the shuttle heat loss in the displacer are given in Fig. 5. 
The energy lost due to internal conduction is negligible in 
the heater and in the cooler and is about 8.5 kW in the 
regenerator, which represents 35% of the total energy loss. 
This is caused by the lengthwise temperature variation, 
which is very significant in the regenerator. The energy lost 
due to dissipation is mainly observed in the regenerator, 
which reaches a maximum of 3.9 kW, with an average of 
935 W. In the heater and in the cooler, it is equal to 26.6 
and 123 W, respectively. The average heat flow value lost 
by the shuttle effect is about 3.1 kW; it represents 13% of 
the total energy loss. 

The energy lost due to external conduction in the 
regenerator is 27 kW, which represents 47% of the total 
losses, as shown in Fig. 6. It is very significant and depends 
mainly on the regenerator efficiency. The energy lost due to 
irreversibility in the compression and expansion spaces is 
very low [18]. 

5. Performance optimization of Stirling engines 

The energy losses are mainly located in the regenerator. 
They are primarily due to the losses by external and 
internal conduction and pressure drop through the heat 
exchangers. The energy lost due to shuttle effect in the 
exchanger piston is also significant; it is about 13%. The 
other losses are very small [17]. 

The reduction of these losses improves the engine 
performance. Such losses depend mainly on the matrix 
conductivity of the regenerator, its porosity, the inlet 
temperature variation, the working gas mass flow rate, the 


regenerator volume, and the geometrical characteristics of 
the displacer. 

To investigate the influence of these parameters on the 
prototype performance, we have changed, each time, the 
studied parameter in the model and have kept others 
unchanged and equal to the prototype parameters. 

5.1. Effect of the regenerator matrix conductivity 
and heat capacity 

The performance of the engine depends on the con¬ 
ductivity and heat capacity of material constituting the 
regenerator matrix. Fig. 7 shows that an increase of matrix 
regenerator thermal conductivity leads to a reduction of 


x 10 S 



Table 2 

Comparison of various model results 


Numerical model 


Heat (J/cycle) Indicated power output Thermal efficiency (%) 


(W) (J/cycle) 


Adiabatic model 

Urielli and Berchowitz [14] adiabatic model 
Dynamic model without losses 
Urielli and Berchowitz [14] quasi-steady flow 
Dynamic model with loss dissipation = (Ml) 

Urielli and Berchowitz [14] quasi-steady flow (pressure drop included) 
(Ml) + Internal conduction loss = (M2) 

M2 + External conduction loss = (M3) 

M3 + Shuttle loss = (M4) 

M4 + losses by gas spring hysteresis = dynamic model with losses 
Experiment 


327 

8286.7 

198.62 

62.06 


8300 


62.5 

314 

7109.3 

170.4 

54.96 


7400 


53.1 

291 

6372.4 

152.47 

53.3 


6700 


52.5 

294 

6355.2 

152.32 

52.64 

314 

6061 

145.27 

46.94 

352 

5886.1 

141 

40.66 

262 

4273 

99.5 

38.49 

— 

3958 

— 

35 


Urieli and Berchowitz model and experimental results are shown in italics. 
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Fig. 5. Lost heat flow in the engine. 
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performance due to the increase of internal conduction 
losses in the regenerator [18]. Fig. 8 shows that the engine 
performance improves when the heat capacity of the 
regenerator matrix rises. 

The matrix of the regenerator can be made from several 
materials. The performance of the engine depends on the 
matrix material. To increase heat exchange of the 
regenerator and to reduce internal losses by conductivity, 
a material with high heat capacity and low conductivity 
must be chosen. Stainless steel and ordinary steel are the 
most suitable materials to make the regenerator matrix. 


5.2. Effect of regenerator porosity 

Porosity of the regenerator is an important parameter 
for engine performance. It affects the hydraulic diameter, 
dead volume, velocity of the gas, regenerator heat transfer 
area, and regenerator effectiveness; thus, it affects the 
losses by external and internal conduction and the 
dissipation by pressure drop [18]. 

The engine performance decreases when porosity in¬ 
creases due to an increase in the external conduction losses 
and a reduction of the exchanged energy between the gas 
and the regenerator ( Q T ), as shown in Fig. 9. 

The performance decreases when porosity increases, but 
porosity has a low limit from which the model does not 
converge, due to the non-satisfaction of the boundary 
conditions. For the prototype studied, the calculated 
optimal porosity is 65.5%, as shown in Table 3. 

5.3. Effect of regenerator temperature gradient: 

(Tf- r —T r _h) 


Although engine losses increase when the temperature 
gradient of the regenerator rises [18], the performance of 
the engine also increases, as shown in Fig. 10. In this event, 
this is due to the increase of energy exchanged between the 
matrix and the working fluid of the regenerator. 

5.4. Effect of regenerator volume 

To vary the regenerator volume, the diameter is fixed 
and the length is varied or conversely. When the regene¬ 
rator diameter is fixed at 0.0226 m, length affects the 
performance. Although the energy exchanged increases, 
engine power and efficiency reach a maximum. When the 
length is equal to 0.01 m, power decreases quickly, as 
shown in Fig. 11 . This can be explained by an increase of 
the dead volume. 

When the regenerator length is constant and L = 
0.022 m, the performance decreases when the regenerator 
diameter increases, as shown in Fig. 12. The dead volume 
and the exchanged energy in the regenerator also decrease. 

5.5. Effect of fluid mass 

An increase of the total mass of gas in the engine leads to 
a rise in the density, mass flow, gas velocity, load, and 
function pressure. Therefore, an increase in the total mass 
of gas in the engine leads to more energy loss by pressure 
drop [18]; however, the engine power increases and the 
efficiency reaches a maximum of about 40% when the mass 
is equal to 0.8 g, as shown in Fig. 13. When the mass 
increases, the decrease of efficiency is due to an increase of 
pressure loss and the limitation of heat exchange capacity 
in the regenerator and the heater. The use of mass of gas 
equal to 1.5 g in the engine leads to an acceptable output 
and a higher power than in the prototype. 
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Fig. 8. Effect of regenerator heat capacity on performance. 


5.6. Effect of expansion volume and exchanger 
piston conductivity 

The expansion volume and the exchanger piston 
conductivity considerably affect the losses due to shuttle 
effect, which represent 13% of the engine total losses. To 
vary the expansion volume, we can maintain the stroke 
constant and vary the piston surface or conversely. 

When the piston stroke is constant and equal to the 
prototype value of 0.046 m, the effect of the piston surface 
on the performance is given in Fig. 14. When the section 
increases, the engine power increases, but the efficiency 
reaches a maximum. 


If the exchanger piston area is equal to 0.0045 m , a 
power higher than 5 kW and an output slightly lower than 
that of the prototype can be reached. When the exchanger 
piston area is constant and equal to the prototype value of 
0.0038 m , the effect of stroke variation on performance is 
given in Fig. 15. When the stroke increases, the engine 
power decreases but the efficiency reaches a maximum. 

The optimal performances are superior to that of the 
prototype. They are obtained when the area and the stroke 
are, respectively, equal to 0.0038 m" and 0.042 m, which 
correspond to a power of 4500 W and an efficiency of 41%. 

The thermal conductivity of the exchanging piston 
affects the engine performances considerably, as shown in 
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Table 3 

Effect of optimal parameters on engine performance 


Result optimization 

Optimal 

value 

Power (W) 

Efficiency 

(%) 

Regenerator 
exchanged energy (J) 

Average pressure 
function (MPa) 

Regenerator porosity 

65.5% 

4554 

43.3 

474.3 

4.72 

Regenerator length (m) 

0.021 

4675 

41 

472.4 

4.79 

Regenerator diameter (m) 

0.024 

4765 

39.9 

475.8 

4.85 

Working fluid mass (g) 

1.15 

4850 

40.1 

480.7 

4.9 

Exchanger piston conductivity W/(mK) 

1.2 

5079 

50.9 

479.9 

4.95 

Exchanger piston area (nr) 

3.86 x 10“ 3 

5183 

50.97 

483.7 

4.92 

Exchanger piston stroke (m) 

0.047 

5106 

51.1 

472.7 

4.95 

Model result before optimization 


4273 

38.5 

448.7 

4.67 

Experiment result of prototype 


3958 

35 
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Fig. 10. Effect of regenerator temperature gradient on performance and exchanged energy. 
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Fig. 16. Weak conductivity reduces the losses by shuttle 
effect and consequently increases the engine power and 
efficiency. 

6. Design optimization of Stirling engine parameters 

Based on the reduction of losses and optimization of 
performance, the optimization consists first in determining 
the optimal value of each parameter when the other 
parameters are equal to the prototype parameters, second 
in gradually replacing the parameters of the prototype by 
their optimal values. The results are presented in Table 3. 


The optimal parameters for the design are found as 
follows. 

For the first line of Table 2, the prototype porosity 
(0.697) is the only parameter replaced by the optimal 
porosity (0.655) in the model. The power and efficiency are 
improved, but the average pressure increases slightly. 
Therefore, the optimal porosity is of the model is 
maintained and the regenerator optimal length (0.021 m) 
is to be searched. By introducing this length into the model, 
the power improves and the efficiency remains acceptable. 
Then, the optimal diameter (0.024 m) is used; the power 
and the exchanged energy in the regenerator increase but 
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Fig. 14. Effect of exchanger piston area on engine performance. 


the efficiency decreases slightly. By replacing the working 
fluid mass by the calculated optimal value (1.15 g) in the 
model, the power and the efficiency increase but the calcu¬ 
lated engine average pressure remains close to the initial 
value. 

The exchanger piston conductivity influences the loss by 
shuttle effect, its reduction leads to a remarkable increase 
in the performances. By introducing optimal values of the 
area and the stroke of the exchanger piston in the model, 
the prototype performance clearly improves. 

7. Conclusion 

The theoretical Stirling cycle has a high theoretical 
efficiency; however, the constructed prototypes have a low 


output because of considerable losses in the regenerator 
and the exchanger piston. This is primarily due to losses by 
external and internal conduction, pressure drops in the 
regenerator, and by shuttle effect in the exchanger piston. 
These losses depend on the geometrical and physical 
parameters of the prototype design. 

An optimization of these parameters has been carried 
out using the GPU-3 engine data, and has led to a 
reduction of losses and to a notable improvement in the 
engine performance. We first applied the parameters of 
this prototype on the developed model; the results were 
very close to the experimental data. Then, we studied 
the influence of each geometrical and physical parameter 
on the engine performance and the exchange energy of 
the regenerator. The reduction of matrix porosity and 
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Fig. 16. Effect of displacer thermal conductivity on performance. 


conductivity of the regenerator increases the performance. 
A rise of the total gas mass leads to an increase of the 
engine power and working pressure. However, the effi¬ 
ciency reaches a maximum. When the displacer section 
increases and the piston stroke decreases, the engine power 
increases, and the efficiency reaches a maximum. A low 
conductivity of the exchanger piston reduces the losses by 
shuttle effect and consequently increases the engine power 
and efficiency. 

Finally, we optimized these parameters gradually by 
introducing them into the model and seeking each time the 
optimum value. Although the actual performance of the 
used prototype is relatively high, results have improved 
considerably. The efficiency increased from 39% to 51%, 


the power rose approximately 20%, and the average 
pressure slightly increased. 

Applying the proposed approach to prototypes with low 
performance or to newly designed engines will lead to the 
determination of their optimal design parameters and 
consequently to a higher performance. 
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